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I. INTRODUCTION 



In the past years, much effort of physicists has been devoted to the understanding of non- 
equihbrium dynamic processes. Phase ordering dynamics, spin glass dynamics, structural 
glass dynamics and interface growth etc are important examples. Since the pioneer work by 
Janssen et al the universal dynamic scaling form in critical dynamics has been explored 
up to the macroscopic short-time regime [l - 10|, when the system is still far from equilibrium. 
Although the spatial correlation length is still short in the beginning of the time evolution, 
the short-time dynamic scaling form is induced by the divergent correlating time around a 
continuous phase transition. Based on the short-time dynamic scaling form, new methods 
for the determination of both dynamic and static critical exponents as well as the critical 
temperature have been developed 0, 0, 1?-!!]. Since the measurements are carried out in 
the short-time regime, one does not suffer from critical slowing down. 

In understanding the dynamic scaling form far from equilibrium, we should keep in mind 
that it holds after a time scale tmic, which is sufficiently long in the microscopic sense, but 
still short in the macroscopic sense. More importantly, the macroscopic initial condition 
should be taken into account in the dynamic scaling form Il2| . For the dynamic 



relaxation starting from an ordered state, i.e., a state with an initial magnetization mo = 1, 
for example, the magnetization decays by a power law , [l2|. If ruQ is smaller but 
close to 1, there emerge corrections to scaling. For the dynamic relaxation starting from 
a random state, i.e., a state without spatial correlations and with a small mo, however, 
the magnetization does not decay, and rather shows an initial increase in the macroscopic 
short-time regime. An independent critical exponent Xq must be introduced to describe the 
scaling dimension of the initial magnetization , [lO|. If mo = 0, the magnetization 



naturally remains zero during the dynamic evolution, but Xq is still needed to describe the 
auto-correlation function etc. This critical exponent also explains the power-law decay of 
the remanent magnetization in spin glasses {2, [isl [l^. On the other hand, the short-time 
dynamic scaling form is universal, in the sense that it does not depend on the microscopic 
details of the dynamic system, such as the lattice types, interactions, and updating schemes 
etc. Up to now, the dynamic relaxation with the ordered and random initial states has been 
systematically investigated. 

Recent progress in the non-equilibrium critical dynamics and its applications includes, for 



2 



example, theoretical calculations and numerical simulations of the XY models and Josephson 



junction^ arrays 
tems 
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15l-ll8|. magnets with quenched disorder |19l-|23|. and various critical sys- 



27l |. Dynamic reweighting methods have been proposed 
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and the dynamic 



29|-|3l|. Recently, the 



approach to the weak first-order phase transitions is also attractive 
ageing phenomenon around a continuous phase transition has been also intensively studied 
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. In this case, the dynamic scaling form for the ageing phenomenon is induced by 
the long-range time correlation, different from that induced by meta-stable states in glassy 
systems below the transition temperature Tc. 

On the other hand, in the past years many activities have been devoted to the domain-wall 
dynamics 
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45(1 . For magnetic materials, for example, a domain wall separates domains 



with different spin orientations. Microscopically, the domain wall may move and create 
bubbles, and macroscopically, it may propagate and roughen. At the zero temperature, there 

3- 



occurs a pinning-depinning phase transition induced by quenched randomness [39|, |46|, 
For a magnetic system with weak disorder, the domain wall does not propagate unless the 
external magnetic field h exceeds a threshold he. At the critical field he, a roughening 
phenomenon is also observed [48]. When a periodic external field h{t) = hn cosjut) is 



applied, the second-order phase transition is softened to a hysteresis loop |39l-l41l. |43| . Most 
these works concerning the domain-wall motion concentrate on the stationary state at the 
zero or very low temperatures and in response to the external magnetic field h{t). 

In this paper, we systematically investigate the dynamic relaxation with a domain wall 
at the critical temperature. For simplicity, we assume that no external field is applied. 
Macroscopically, therefore, the domain wall does not propagate. We should only keep in 
mind that different from the case at the zero temperature, here the bulk also evolves in 
time. To be specific, we consider the dynamic relaxation starting from a semi-ordered initial 
state. For the Ising model, for example, the semi-ordered initial state consists of two fully- 
ordered domains with opposite spin orientations. As time evolves, the domain wall roughens, 
and looks like a growing interface. In this paper, we call it the domain interface. Such a 
domain- wall dynamics is theoretically and practically important. Theoretically, it is very 
interesting to investigate the short-time dynamic scaling behavior starting from the semi- 
ordered state, in comparison with that starting from the ordered or random state. It extends 
the study of the domain-wall motion at the zero or very low temperatures to the critical 
temperature, and especially explores the dynamic behavior far from equilibrium. In this 
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paper, we intend to clarify first the dynamic scaling behavior of model A 



Then the 



dynamic theory may be generalized to model B. Along this direction, one mig ht find the 



way to study relevant dynamic processes of driven diffusive lattice gases |24l. |50|. 

Furthermore, the non-equilibrium critical dynamics around a surface is also an important 
topic 
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54| . For the dynamic relaxation starting from the random state, the dynamic 
evolution of the magnetization at surface is controlled by both the scaling dimension xo 
of the global initial magnetization and the static exponent /3i of the surface magnetization 
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52| . For the dynamic relaxation starting from the ordered state, it is expected that /3i is 
sufficient to describe the dynamic evolution of the magnetization at surface. For the dynamic 
relaxation without a surface but starting from the semi-ordered state, it looks somewhat like 
that there exists a fictitious surface. The dynamic evolution of the magnetization inside the 
domain interface is governed by an exponent Pi. But this /3i does not correspond to a static 
exponent in equilibrium, and it is induced by the semi-ordered initial state. Therefore, an 
additional purpose of this paper is to compare the dynamic relaxation starting from the 
semi-ordered state with that starting from the ordered state but around a surface. 

In Ref. 55I, brief results on the dynamic relaxation of the magnetization have been 
reported for the two-dimensional Ising model. This paper aims at a comprehensive study 
of the topic, and explores especially the dynamic scaling behavior of the Binder cumulant 
(or susceptibility), height function and roughness function. Furthermore, Monte Carlo sim- 
ulations are performed also for the three-dimensional Ising model, to study the dimension 
dependence of the scaling functions and critical exponents as well as the corrections to scal- 
ing. In Sec. II, the models and scaling analysis are described, and in Sec. Ill, the numerical 
results are presented. Finally, Sec. IV includes the conclusions. 

II. MODEL AND SCALING ANALYSIS 
A. Model 

The d- dimensional Ising model is the simplest model for magnetic materials, exhibiting 
a second-order phase transition. The Hamiltonian is written as 

-^h = ky: s.s, , (1) 
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where Si = ±1 is an Ising spin at site i of a square or cubic lattice, the sum is over the nearest 
neighbors, and T is the temperature. In this paper, we set the temperature at its critical 
value Tc] or in other words, we set K at its critical value Kc- The Hamiltonian of the Ising 
model itself does not include an intrinsic dynamics. For example, Monte Carlo algorithms 
may be introduced to simulate the dynamic evolution of the system. It is generally believed 
that the Monte Carlo dynamics is in the same universality class of the Langevin equation. 

Let us consider a kind of dynamic relaxation processes at the critical temperature. After 
a macroscopic initial state at very low temperatures is prepared, the dynamic system is 
suddenly quenched to the critical temperature, and then released to the dynamic evolution 
of model A j^, I49I. For the dynamics of model A, the order parameter and other relevant 
physical quantities are not conserved during the dynamic evolution. In Monte Carlo simu- 
lations, it can be simply realized with a standard one-spin flip. In this paper, the heat-bath 
algorithm is always used in the dynamic Monte Carlo simulations. Selecting a single spin 
Si, we flip it with the transition rate. 



where j{i) labels the nearest neighbors of the site i, and c is the normalization constant. In 
fact, other Monte Carlo algorithms, such as Metropolis algorithms, Monte Carlo algorithms 
with a multi-spin flip and rejection-free Monte Carlo algorithms etc, yield the same results. 
The condition is that the algorithms should be local, i.e., only spins in a local region are 
flipped in a single flip. 

With Monte Carlo simulations, we first study the critical relaxation starting from a semi- 
ordered state, taking the two-dimensional (2D) and three-dimensional (3D) Ising model as 
examples. The Ising model is defined on a rectangular lattice 2L x L in two dimensions and 
2L X in three dimensions, with a linear size 2L in the x direction and L in the other 
directions. Periodic boundary conditions are used in all the directions. The semi-ordered 
state is such a state, that spins are positive on the sublattice L'^ {d = 2 or 3) at the right 
side and negative on the sublattice L'^ at the left side. For convenience, we set the x-axis 
such that the domain wall between the positive and negative spins is located at x = 0. So 
the X coordinate of a lattice site is a half-integer. 

After preparing the semi-ordered initial state, we update the spins with the heat-bath 
algorithm at the critical temperature Tc. Since no external magnetic field is added, macro- 




exp{K J2j{i) Sj) + exp{-K J2j{i) S. 



(2) 
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scopically the domain wall does not move. As time evolves, however, the domain wall 
fluctuates and creates bubbles. As a result, the domain wall becomes thicker and thicker, 
and a dynamic roughening process occurs. Therefore, we call it a domain interface. In Fig.[T], 
the dynamic evolution of the spin configuration around the domain wall is illustrated. Some- 
what different from a standard growing interface, here the bulk evolves in time. In analyzing 
the dynamic properties of the domain interface, this must be kept in mind. 

For comparison, we also perform Monte Carlo simulations of the Ising model with a free 
or disordered surface, but starting from the ordered state. In this case, the lattice is taken to 
be L'^ {d = 2 or 3). For the free surface, a free boundary condition is used in the x direction, 
while periodic boundary conditions are used in other directions. For convenience, we set the 
X-axis such that the free surface locates at x = 1 or L. For the disordered surface, the spin 
Si on the surface couples to a random field hi through the interaction —Hi/kT = KhiSi 
with the random filed hi = ±1. The disordered boundary condition is only implemented in 
the X direction, and periodic conditions are used in other directions. Since the initial state 
is the ordered state, the magnetization decays in time. In particular, the magnetization is 
also x-dependent due to the geometric surface. The dynamic behavior of the magnetization 
around the surface is governed by the surface exponents, while that at bulk is controlled 
by the bulk exponents. The region affected by the geometric surface grows in time, and it 
looks like that the surface becomes thicker and thicker. Phenomenologically, this dynamic 
behavior is similar to that of the domain interface. 

Finally, to expose the dynamic evolution of the bulk, we perform Monte Carlo simulations 
of the Ising model with periodic boundary conditions in all directions, starting from the 
ordered state. The lattice is taken to be L'^ in d dimensions. 

B. Scaling analysis 

We first analyze the dynamic scaling behavior of the domain interface. Due to the semi- 
ordered initial state, the time evolution of the dynamic system is inhomogeneous in the x 
direction. Therefore, we measure the magnetization and its second moment as functions of 
X and t. In two dimensions, for example. 




k 



1,2. 



(3) 
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Here Sxyif) is the spin at the time t on the lattice site {x,y), L is the lattice size, and 
< . . . > represents the statistical average. For convenience, we also use M{t, x) = M^^Ht, x) 
to denote the magnetization. Then we can define a time- dependent Binder cumulant jol. Issl. 

U(tx) = ^^^-1 (4) 

The susceptibility M(2)(t, x) - M{t, x)^ or the Binder cumulant U{t,x) describes the fiuctu- 
ation in the y direction. In three dimensions, we simply use S^yz to denote the spin on the 
lattice site {x,y, z), and similarly define the magnetization and Binder cumulant. 

In order to characterize the growth of the domain interface and its fiuctuation in the x 
direction, we introduce a height function and its second moment in the x direction. 



A; = 1,2. (5) 



lx=l J / 

Here < . . . > represents not only the statistic average but also the average in the y direction. 
As usual, we also use the notation h{t) = h^^^{t). Then the roughness function of the domain 
interface is defined as 

uj\t) = h^^\t)-h{t)h{t). (6) 

Except for the scaling dimension of the magnetization, the height function measures the 
thickness of the domain interface, while the roughness function represents the fiuctuation of 
the domain interface. 

At the critical temperature Tc and in the thermodynamic limit, there are two length 
scales in the dynamic system, i.e., x and the non-equilibrium spatial correlation length ^(t). 
For a finite system, the lattice size L is an additional length scale. In general, one may 
believe that ^(t) universally grows as ^{t) ~ t^^^ in all spatial directions, because of the 
homogeneity of the interactions in the Hamiltonian. Therefore, general scaling arguments 
lead to the scaling form of the magnetization and its second moment 

M(^)(t, X, L) = t-^^''''M^^\t^'^/x, t^^'/L), k = l,2. (7) 

Here f3 and u are the static exponents, and z is the dynamic exponent. On the right side 
of the equation, the overall factors f'^PI'^^ indicates the scaling dimension of M'^^\ and the 
scaling function M*^'^'^(t^/^/x, t^/^/L) represents the scale invariance of the dynamic system. 



In general, the scaling form in Eq. 
after a microscopic time scale tmic 111, 



holds already in the macroscopic short-time regime, 

a. 



For the magnetization, the scahng function M{i}l^ jx^t^l^ jV) is independent of L in the 
thermodynamic hmit L ^ oo. Then the scahng form is simphfied to 

M(t,x) =t-^/""M(t^/7x). (8) 

For the susceptibihty, it is different. For a sufficiently large lattice and in the short-time 
regime, the non-equilibrium spatial correlation length ^(t) is much smaller than the lattice 
size L. Therefore, the spatially correlating terms < S^yj^Sxy^ > — < Sxy^ >< Sxy2 > with 
I2/2 — yi\ > ^{t) can be neglected. In other words, one of the two summations over yi and 
y2 in the susceptibility M^'^^(t,x) — M{t^xY is suppressed. It then leads to the finite-size 
behavior M^'^\t,x) - M{t,xf ~ l/L'^-'^ (rf = 2 or 3) j9|]. Together with Eqs. ^ and 
one may derive the scaling form of the Binder cumulant 

U{t,x) = ^j^U{t'/yx). (9) 

The Binder cumulant is interesting, for only the dynamic exponent z is involved. 

By definition, the height function h{t) is nothing but the average magnetization in the 
positive domain, i.e., h(t) = J2x>o ^(t^^)/ In general, h{t) does not obey a simple power 
law. Its behavior replies on the scaling function M(t^^^ /x). In fact, one may deduce a 
scaling form for h{t) from Eq. ([7]), 

h{t) = r^/^'hit^/yL). (10) 

Different from M{t,x), here one should not ignore the dependence on the lattice size L, for 
the scaling function h{i}-/^ /L) just represents the dynamic effect of the domain interface. 
This is obvious from the definition h{t) = J2x>o^{'tj^)/L- On the other hand, one should 
also note that the height function here is scaled as the magnetization, not a spatial length. 

Similar to Eq. ( ITOl) . one may also assume that the scaling form for the roughness func- 
tion is u}'^{t) = t~'^^/'^/^F{t}/^ / L). For later convenience, we separate a factor t^l^ jL from 
F{i}l^ I L\ and rewrite the scaling form as 

A\-2filv)lz _ 

uj\t) = — - — u^t'/yL). (11) 

In general, uj'^{t) does not exhibit a power-law behavior. This is different from a standard 
growing interface. The reason is that here uj'^{t) includes fluctuations from the domain 
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interface and the bulk. In fact, we will show in the next section that the scaling function 
up-ii}!^ jV) describes the fluctuation induced by the domain interface. 

The scaling forms in Eqs. (I7|) — f llip can be also applied to the dynamic relaxation with 
a free or disordered surface, but starting from the ordered state. One should only keep in 
mind that /3, v and z are the critical exponents at bulk. The critical exponents at surface 
should be deduced from the scaling functions in Eqs. (17|1 — (ITT]) . In this case, the lattice is 
taken to be L'^. The inhomogeneity in the x direction is induced by the surface. 

The purpose of this paper is to investigate whether the scaling forms in Eqs. (|7|) — (ITT!) do 
hold in the dynamic relaxation with the domain interface and with the free or disordered 
surface. With Monte Carlo simulations, we study characteristics of the scaling functions, 
and extract corresponding critical exponents. Dynamic systems with the domain interface 
and with the free or disordered surface share some common features, although they are 
intrinsically different. The domain interface is induced by the geometric structure of the 
semi-ordered initial state, while the dynamic relaxation with the free or disordered surface 
is controlled by the geometric surface which remains even in equilibrium. 

It is important that the height function h{t) and roughness function in Eqs. ([5]) and 
([6]) include the dynamic evolution of the bulk. Therefore, their behaviors deviate from those 
of a standard growing interface. To obtain the dynamic features of a pure growing interface 
such as the power-law behavior and the roughness exponent etc, we need to subtract the 
contribution of the bulk. Therefore, we finally perform Monte Carlo simulations of the Ising 
model on a lattice L'^ with periodic boundary conditions in all directions, and starting from 
the ordered state. In this case, the dynamic system is homogeneous in all directions. The 
height function hh{t) and the roughness function Co'^(t) are just the line magnetization and 
its susceptibility in the x or y direction. The scaling functions h{t}/^ /L) and uj'^it^/^ / L) in 
Eqs. (fTOj) and (fTTj) are constants. In other words, hi,{t) and show a power- law behavior. 

Then we may redefine the pure height function and roughness function for the dynamic 
relaxation of the domain interface or around the surface by subtracting the contribution 
from the bulk 



We may expect that Dh{t, L) and Du^{t, L) exhibit a power-law behavior as in the case of 



Dh{t,L) = hb{t) - hit), 



(12) 



Duj''{t,L)=uj\t)-ujl{t). 



(13) 
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a standard growing interface. Here we should note that we define Dh(t,L) as hb(t) — h{t) 
rather than h[t) — hb{t), for h{t) decays in time faster than hh{t). 



III. MONTE CARLO SIMULATION 



For the 2D Ising model, our main results are presented with L = 512 at = 0.44069, 
and the maximum updating time is tM = 25600. Additional simulations with L = 1024 
and L = 256 are performed, to investigate the finite-size scaling behavior and finite-size 
effect. The total samples for average are 24000. For the 3D Ising model, the main results 
are obtained with L = 128 at Kc = 0.22165, and the maximum updating time is Im = 2560. 
Additional simulations with L = 256 and L = 64 are performed to investigate the finite- 
size scaling behavior and finite-size effect. The total samples for average are 30000. The 
statistical errors are estimated by dividing the total samples into two or three subgroups. If 
the fluctuation in the time direction is comparable with or larger than the statistical error, 
it will be taken into account. 

Theoretically, the scaling forms described in the previous section hold in the macroscopic 
short-time regime, after a microscopic time scale tmic- tmic is not universal, and relies on 
microscopic details of the dynamic systems. In Monte Carlo simulations, for example, t^ic 
is typically tens or hundreds of Monte Carlo time steps [9]. With quenched disorder or 
frustration, tmic could be longer. For the simple Ising model with the nearest neighbor 
interactions, tmic is rather short, about 10 - 20 time steps. Therefore, critical exponents 
are typically obtained in the time intervals [20, 25600] in two dimensions and [10, 2560] in 
three dimensions. From the data collapse of different x and L, one may observe the scaling 
functions in a even larger time window. 



A. Magnetization 

The time evolution of the magnetization of the 2D Ising model starting from the semi- 
ordered state is displayed in Fig. [2J For a sufficiently small s, e.g., x = 255.5 and t < tM = 
25600, M{t,x) approaches the non-linear decay at bulk, M{t,x) ~ 

t-l3/uz jgj_ rjj^g exponent 

(S/uz = 0.0580(3) measured from the slope of the curve is well consistent with /3 = 1/8, 



1^ = 1 and z = 2.16(1) reported in the literature jsj. For a sufficiently large s, e.g., x = 0.5 
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and t > 20, M{t, x) exhibits also a power-law behavior, but decays much faster than at bulk. 
In other words, we catch some features of the scaling function M{s) in Eq. ([8]), 

— const s — )• 

M(.)~{ , (14) 

with s = t^^^ /x. In the limit s — ?■ oo, one may define an interface exponent /3i such that 

M{t, x) ~ t'^'/''' ■ /3i = /3 + /3o- (15) 

Inside the interface, the power-law decay of the magnetization is governed by the interface 
exponent (3i, while outside the interface, it is controlled by the bulk exponent f3. In Fig. [2], one 
measures (ii/vz = 0.518(4), and then calculates (3i = 1.119(9) and f3o/h' = 0.994(9). Similar 
to the exponent xq in the dynamic relaxation starting from the random state [ll, [a] , Po here 
is induced by the semi-ordered initial state. Accounting the error, one may believe Po/iy = 1, 
which suggests that M(t, x) is an analytic function of x. This result is also supported by 
the simulations of the 3D Ising model. Since /3i is much bigger than /3, the magnetization 
inside the domain interface decays much faster than that at bulk. This phenomenon is 
understandable, for the dynamic evolution of the spins in the positive domain is strongly 
affected by those in the negative domain, and vice versa. 

To fully confirm the scaling form in Eq. ([8]), for example, we fix x' = 1.5, and change 
the time scale t of another x to (x'/x)^ t, and the scale of M(t,x) to {x' / x)~^/'^ M{t,x). As 
shown in Fig. [2], all data of different x nicely collapse to the curve of x' = 1.5. This validates 
Eq. ([8]). Alternatively, we may plot t^/^'M^t function of s = t^/^ /x. According to 

Eq. ([8]), all data of different x should collapse onto the master curve M{s). This is shown 
in Fig. [3l Clearly, M{s) — )■ const when s — 0, while M{s) — )■ s~^o/^ when s — ?■ oo. 

For comparison, the time evolution of the magnetization of the 2D Ising model with a 
free surface but starting from the ordered state is shown in Fig. |H For a sufficiently small 
s, e.g., X = 256 and t < tM = 25600, M{t,x) approaches also the non-linear decay at 
bulk, M{t,x) ~ t~^/'^^ with l3/uz = 0.0579(4). For a sufficiently large s, e.g., x = 1 and 
t > 20, M{t,x) exhibits also a power-law behavior. Assuming again the scaling ansatzes 
in Eqs. f[T^ and f[T^ . the measurement of the slope yields Pi/vz = 0.231(1). Then one 
calculates (3i = 0.499(2). It is in good agreement with the surface exponent /3s = 1/2 for 
the free surface [56|. Now, the exponent /^o/'^ is estimated to be 0.374. Therefore, M(t,x) 
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is not an analytic function of x, when x approaches the free surface. This is very different 
from the domain interface. Finally, the data collapse according to Eq. (jH]) is also shown in 
Fig. m and it fully confirms the scaling form. 

As in Figs. [2] and IH similar analysis can be carried out for the magnetization of the 
2D Ising model with a disordered surface. For a small s, one measures (3/i>z = 0.0578(7). 
For a large s, careful analysis shows that the power-law behavior is not perfect jsSj. In 
the equilibrium state, one may show that the surface exponent /3s of the disordered surface 
remains 1/2, but with a logarithmic correction to scaling 56|. Therefore, we fit the time- 
dependent magnetization at a; = 0.5 with a logarithmic correction to scaling, i.e., M{t) = 
Cit-"V(1 + C2ln(^))^/^ and derive ai = 0.231, consistent with Pi/uz = 0.231(1) for the free 
surface. If one fixes C2 = 0, it yields ai = 0.272, significantly different from 0.231(1). 

In Fig. [3|, the scaling function M(s) with s = t^/^ /x is plotted for the domain interface, 
free surface and disordered surface. We clearly observe the characteristic of the scaling 
function in Eq. ( !T4l) . and measure the exponent /3o = 0.998(5) for the domain interface, 
and 0.372(6) for the free surface. Due to the logarithmic correction to scaling, M{s) of the 
disordered surface decays faster than that of the free surface at the large s regime. 

We emphasize that in the case of the free surface or disordered surface, the exponent 
/3i = does describe the critical behavior of the magnetization at the surface in equilibrium. 
Around the free surface, for example, M (r) ~ (— t)^^ with r being the reduced temperature. 
It is important that /3i is induced by the geometric surface which remains forever. In the 
case of the domain interface, however, /3i is induced by the geometric structure of the semi- 
ordered initial state. When the dynamic system reaches its equilibrium state, the influence of 
the initial state disappears and the critical behavior of the magnetization is governed by the 
bulk exponent /3 everywhere. We should keep in mind, however, that exactly at the critical 
temperature (i.e.,r = 0) and in the thermodynamic limit, the dynamic system never reaches 
its equilibrium state in a finite time due to the divergent correlating time. According to 
Eq. (IHl), therefore, the domain interface and the free or disordered surface behave similarly. 

For the 3D Ising model, the static and dynamic exponents at bulk are known to be 
/3 = 0.327(1), z/ = 0.630(2) and z = 2.04(1) [57|. For the free surface, the surface exponent 
is (3s = 0.795(10) 581]. Following the procedure for the 2D Ising model, we have analyzed 
the scaling behavior of the dynamic relaxation with the domain interface, free surface and 
disordered surface. Especial attention is drawn to the critical exponent Po/iy. 
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Let us first consider the domain interface. For a small s, the magnetization shows the 
power-law behavior at bulk, M(t, x) ~ t~^^'^^. The critical exponent is estimated to be 
/S/uz = 0.253(5), well consistent with the value 0.253(1) at bulk js^. For a large s, e.g., 
X = 0.5 and t > 10, the magnetization exhibits the power-law behavior M(t, x) ~ f-Pi/^z 
with /?! = /3 + /3o in Eq. ( ITB]) . From the slope of the curve, one obtains Pi/uz = 0.744(2). 
Then one calculates the critical exponent (3o/u = (0.744 - 0.253) x 2.04 = 1.002(4). For 
the 2D Ising model, (3o/h' = 0.998(5). These two measurements of (3o/i^ strongly suggest 
Po/'^ = 1, and it is dimension- independent. Therefore M{t,x) is an analytic function of x. 
In Fig. 13 the scaling functions M{s) of the magnetization for the 3D Ising model is plotted. 
Data collapse for different x is observed. From the slope of the curve in the large s regime, 
one measures Pq/u = 1.001(6). 

Similar analysis is applied to the magnetization of the 3D Ising model with the free and 
disordered surfaces, and the scaling function is also shown in Fig. [51 Different from the 
case of the 2D Ising model, the large-s tails of the scaling function M{s) for the free and 
disordered surfaces look parallel each other. In the inset of the figure, the magnetization at 
X = 1.0 is displayed for both the free and disordered surfaces. The slope is Pi/vz = 0.623(5) 
for the free surface, and 0.632(2) for the disordered surface. The difference is only one or two 
per cent, and the correction to scaling is rather small. For the free surface, one estimates 
Pi = 0.623(5) X uz = 0.801(6). Alternatively, one may also measure Po/i^ = 0.747(6) from 
the scaling function obtained with different x, and then calculates /3i = /3 + /3o = 0.798(4). 



These values of /3i are well consistent with the surface exponent /3s = 0.795(10) (58|. 



B. Binder Cumulant 



In Fig. El the time evolution of the Binder cumulant is displayed for the 2D Ising model 
starting from the semi-ordered state. For a sufficiently small s, e.g., x = 255.5 and t < tM = 
25600, the Binder cumulant exhibits the power-law behavior at bulk, U{t, x) ~ From 
the slope, one measures {d — l)/z = 0.468(4), and it is consistent with {d — l)/z = 0.463(3) 
calculated from z = 2.16(1). For a sufficiently large s, e.g., x = 0.5 and t > 20, the Binder 
cumulant grows also by a power law, but much faster than that at the large x. Then we 
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extract the characteristic of the scahng function U{s), 



const 



s ^ 



U{s) ~ { 



(16) 



s' 



s — )■ oo 



In the hmit s — )■ oo, one may derive from Eqs. ([9]) and ( fT6l) 



(17) 



From the curve of x = 0.5 in Fig. [6l one estimates {d — l + dQ)/z = 1.390(8). Taking z = 2.16 
as input, one calculates cIq = 1.390 x 2.16 — (rf — 1) = 2.00(2), very close to 2. 

In Fig. [71 the Binder cumulant is plotted for the 3D Ising model starting from the semi- 
ordered state. For a small s, one observes the power-law behavior at bulk, U{t, x) ~ 
From the slope of the curve, one obtains 0.995(12), and then estimates z = 2.01(2), consistent 
with z = 2.04(1) from the literature [57]. For a large s, e.g., x = 0.5 and t > 10, one estimates 
{d—l+dQ)/z = 1.963(10) from the power-law behavior in Eq. ( fTTI) . then derives do = 2.01(2). 
Again it is close to 2. To further verify the scaling form in Eq. (|9]), we fix x' = 1.5, and 
change the time scale t of another x to {x' /xY t, and the scale of U {t, x) to [x' / xY^^U {t, x). 
As shown in Fig. [6]and[7l all curves of different x nicely collapse to the curve of x' = 1.5. 
This fully confirms the scaling form in Eq. (Q. 

To reveal the lattice-size dependence of the Binder cumulant in Eq. (I9l), we fix x = 0.5, 
and plot U{t,L) = U{t,x = 0.5) as a function of t for different L in Fig. [HI Obviously, all 
curves of different L and in two- and three-dimensions are parallel each other, and exhibit 
the power-law behavior in Eq. f[T71) . We then fix a lattice size L', and change the scale of 
U{t,L) of another L to U {t, L){L' / L)~^'^~'^\ Data collapse is clearly observed for both the 
2D and 3D Ising models. 

In Fig. [9], the scaling function U{s) with s = t^^^/x is plotted for the 2D Ising model 
with the domain interface and free surface. Data of different x collapse clearly onto their 
master curves. For the domain interface, the asymptotic behavior of U{s) in Eq. f[T6|) is 
exposed. The exponent rfg is measured to be 2.00(2), the same as that extracted from the 
single curve of x = 0.5 in Fig. [6l For the free surface, U{s) const is also observed in 
the limit s — 0. In the large s regime, however, it does not exhibit a power-law behavior. 
Instead, it increases by a logarithmic law, U (s) = aQ + ai ln(s). In other words, the exponent 
do of the free surface is effectively but with a logarithmic correction. This result indicates 
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that the spatial fluctuation of the domain interface grows in time much faster than that of 
the free surface. 

In three dimensions, rf — 1 — 2/3i/z/of the free surface is negative. Starting from an ordered 
state, the susceptibihty decays in time. Therefore, one suffers from large fluctuations, and 
it is difficult to address the dynamic behavior of the susceptibility or Binder cumulant. For 
the disordered surface, the situation is even more complicated. Since our paper is already 
lengthy, we will not go into the details here. 



C. Height function and roughness function 

In the preceding two subsections, we have analyzed the temporal and spatial structures 
of the magnetization M{t,x) and Binder cumulant U{t,x). Up to now, however, we have 
not yet touched how the interface grows and fluctuates in the x direction. For this purpose, 
we have introduced the height function h{t) and the roughness function in Eqs. ([5]) 

and For a standard growing interface, the time evolution of the height function h{t) 
may be not so important, but the roughness function Co'^(t) exhibits a power-law behavior 
governed by the roughness exponent. 

In Fig. [ini the height function h{t) is plotted for the 2D Ising model. One flnds a power- 
law behavior for the dynamic relaxation of the bulk, i.e., with periodic boundary conditions 
in all directions. The slope of the curve is 0.0576(3), consistent with P/uz = 0.0579(3) from 
the literature 9(]. For the domain interface and free surface, the height function decreases 
faster than a power law. Actually, the curves can be fltted by a double power law, e.g., 
h(t) = cot"° — cit"^ Although this four-parameter flt could not produce very accurate values 
of ao and ai, it leads us to introduce the pure height function Dh{t,L) in Eq. f|T2|) . The 
conjecture is that the term Cit"^ represents the pure interface, and Cot"" is the magnetization 
of the bulk. In Fig. [TUl we do observe a power-law behavior for the pure height function 
Dh{t,L). The slope of the curves is estimated to be 0.407(2). In comparison with that for 
a standard growing interface, this power-law behavior is special for the domain interface. 

In Fig. [m the roughness function is plotted for the 2D Ising model. In the case of the 
bulk, u)'^{t) is expected to obey a power law, although there exist corrections to scaling. 
Anyway, one may roughly estimate the exponent (1 — 2/3/z/)/z to be 0.345(14), consistent 
with z = 2.16(1). Due to corrections to scaling, the dynamic behavior of cj^(t) looks unclear 
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for the domain interface and free surface. However, the pure roughness function Du'^{t,L) 
in Eq. f|T3|) obviously obeys a power-law behavior for both the domain interface and free 
surface. It seems that cj^(t) and uj'^{t) have the same correction to scaling. Neglecting the 
corrections to scaling, one may assume 

~ c without interface , , 

u;\u) = { . (18) 

c + u with interface 



Then one derives 



4^{l-2l3/iy+d^)/z 

Du^%L)= . (19) 



In Fig. [m the curves of Duj^{t,L) for the domain interface and free surface are parallel 
each other. From the slopes of the curves one measures (1 — 2/3/z/ + d^)/z = 0.808(2). 
Then one calculates the exponent d^^ = 0.995(4). The fluctuations of the domain interface 
and free surface grow faster in time than that of the bulk. More importantly, the exponent 
(1 — 2(3/ u + d^)/z is independent of the interface exponent or surface exponent f3i. 

One may also verify the lattice-size dependence, Duj'^{t,L) ~ in Eq. ( fT9|) . In 

Fig. [121 Du'^{t,L) is plotted for different lattice sizes. Obviously, all the curves are parallel 
each other. Then we fix a lattice size, e.g., L' = 256, and change the scale of Du)'^{t,L) 
of another L to Duj'^{t, L)(L' / L)~^^'^'^"\ Taking rf^ = 1 as input, data collapse is clearly 
observed. 

For a standard growing interface, the roughness function grows by a;^(t) ~ t^"/^, and a 
is the so-called roughness exponent. According to Eq. ( fT9l) . the roughness exponent for the 
domain interface and free surface is a = (1 — 2/3/z/ + d^^)/2 = 0.872(8). On the other hand, 
from the dimension counting one may expect Dh{t, L) ~ t"/^ for the pure height function. 
From the measurement a/ z = 0.407(2) in Fig. [TOl one calculates a = 0.879(6). These two 
measurements of the exponent a are in good agreement with each other. In the scaling 
analysis of Duj'^{t,L), —fi/v represents the scaling dimension of the magnetization. One 
may remove it, e.g., by dividing Du)'^{t,L) by M{tY of the bulk. Then the real roughness 
exponent is (1 + di^)/2. Since rf^ = 1, the roughness exponent is just 1. This conclusion 
holds also for the Ising model in three dimensions. But the measurements of the exponents 
are somewhat complicated in this case, for 1 — 2(5/ v = —0.038 is around zero. 
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IV. CONCLUSION 



In summary, we have investigated the non-equihbrium critical dynamics with a domain 
interface, a free surface and a disordered surface, taking the two- and three-dimensional Ising 
models as examples. The dynamic scaling behavior is revealed, and a dynamic roughening 
process is observed. Critical exponents characterizing the magnetization. Binder cumulant, 
height function and roughness function are extracted, and the results are summarized in 
Table [11 

i) For the domain interface, (^o/v for the magnetization in Eq. f|T^ takes values close to 
1 for both the two- and three-dimensional Ising models. It indicates that the magnetization 
M(t, x) is an analytic function of x. Especially, M(t, x) inside the domain interface decays 
much faster in time than that at hulk, for the interface exponent /3i = /3 + /3o is much bigger 
than the bulk exponent /3. For the free surface, the values of /3i = fig are in agreement with 
the measurements in equilibrium. For the disordered surface, /3i = /3* takes the same value 
as that of the free surface, but with a logarithmic correction to scaling in two dimensions. 

ii) For the domain interface, the exponent for the Binder cumulant in Eq. flTB]) takes 
values close to 2 in two and three dimensions. For the free surface, d^ = dg is effectively 
in two dimensions, but with a logarithmic correction to scaling. These results indicate that 
the fluctuation in the y direction inside the domain interface is stronger than that around 
the free surface. In fact, one can derive from Eqs. (IT^ and (fT6|) that inside the domain 
interface, the susceptibility behaves as M^'^\t, x) — M(t, same as that 
at bulk. Around the free surface, the susceptibility is M(t,x)''^^ — M(t, 

different from that at bulk. 

iii) For both the domain interface and free surface, the roughness function in Eq. ffTTj) 
does not obey a power law, for it includes the fluctuation of the bulk and domain interface. 
After subtracting the contribution of the bulk, the pure roughness function in Eq. f|T3l) does 
exhibit a power-law behavior in Eq. ( JT9l) . and the roughness exponent is identified to be 
a = (1 + (i^)/2. Interestingly, the exponent di_j takes values close to 1 for both the domain 
interface and free surface, and also independently of the spatial dimension. In other words, 
the fluctuation of the interface in the x direction is independent of the exponent (3i. 

Theoretically, above results need further understanding. For example, it is a challenge to 
derive the dynamic scaling forms with renormalization group methods. It is also important 
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to investigate how the quenched disorder may affect the domain-wall motion at the critical 
temperature. The techniques used in this paper may be also applied to similar dynamic 
systems. 
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0.630(2) 


[57J 




z 


2.16(1) [gj 


2.04(1) 


[57] 



TABLE I: In the upper sector, the exponents Po/u, j3i and do are for the domain interface, Ps 
and dg are for the free surface, and /?* is for the disordered surface. The exponent is for both 
the domain interface and free surface. In the measurements of /3* and dg for the 2D Ising model, 
logarithmic corrections to scaling are taken into account. In the lower sector, the static exponents 
P and V, the dynamic exponent z, and the surface exponent Ps are taken from the literatures. 
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FIG. 1: Dynamic relaxation from a semi-ordered state is simulated for the 2D Ising model at the 
critical temperature. The spin configuration of the domain interface is shown in a spatial widow 
[-40,40] at the time t = 0,10,100,1000 (from left to right). Black points denote Si = -1 and 
white points denote Si = 1. The lattice size L = 256 is used in the Monte Carlo simulations. 
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FIG. 2: The magnetization of the 2D Ising model starting from the semi-ordered state is plot- 
ted with solid lines on a double-log scale. Dashed lines show the power-law fits. According 
to Eq. ([8]), data collapse for different x is demonstrated at a fixed x' = 1.5. Solid circles, 
solid squares, solid triangles, pluses, open circles, open squares and open triangles correspond 
to X = 0.5, 1.5, 3.5, 7.5, 15.5, 31.5 and 63.5 respectively. From Ref. [55], this figure is reproduced by 
permission of Europhys. Lett.. 
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FIG. 3: The scaling functions M{s) with s = t^/^/x is plotted on a double-log scale, for the 2D 
Ising model with the free surface, disordered surface and domain interface (from above). Data 



collapse for different x is observed. Dashed lines show the power-law fits. From Ref. 
figure is reproduced by permission of Europhys. Lett.. 
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FIG. 4: The magnetization of the 2D Ising with a free surface and starting from an ordered state, 
is plotted with solid lines on a double-log scale. Dashed lines show the power-law fits. According 
to Eq. dSD, data collapse for different x is demonstrated at a fixed x' = 2. Solid circles, solid 
triangles, pluses, open circles, open squares and open triangles correspond to x = 1, 4, 8, 16, 32 and 
64 respectively. 
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FIG. 5: The scaling function M{s) with s = t^^^ /x is plotted on a double-log scale, for the 3D Ising 
model with the free surface, disordered surface and domain interface (from above). Data collapse 
for different x is observed. Dashed lines show the power-law fits. No logarithmic correction is 
detected for the disordered surface. In the inset, the magnetization at x = 1.0 for the disordered 
surface and free surface are shown. 




FIG. 6: The Binder cumulant of the 2D Ising model with the domain interface is plotted with 
solid lines on a double-log scale. Dashed lines show the power-law fits. According to Eq. ([9]), data 
collapse for different x is demonstrated at a fixed x' = 1.5. Solid squares, solid triangles, pluses, 
open circles, open squares and open triangles correspond to x = 0.5,3.5,7.5,15.5,31.5 and 63.5 
respectively. 
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FIG. 7: The Binder cumulant of the 3D Ising model with the domain interface is plotted with 
solid lines on a double-log scale . Dashed lines show the power-law fits. According to Eq. data 
collapse for different x is demonstrated at a fixed x' = 1.5. Solid circles, solid triangles, pluses, open 
circles, open squares and open triangles correspond to x = 0.5, 3.5, 7.5, 15.5, 31.5, 63.5 respectively. 




FIG. 8: The Binder cumulant at x = 0.5 for the Ising model with the domain interface is plotted 
with solid lines on a double-log scale. The lower three solid lines are for the 2D Ising model 
with the lattice size L = 128, 256 and 512. The x— and y— axis are on the bottom and right 
sides. According to Eq. Q, the curves collapse to the curve at L' = 256. Circles and triangles 
correspond to L = 512 and 128 respectively. The upper three solid lines are for the 3D Ising model 
with L = 32, 64 and 128. The x— and y— axis are on the top and left sides. Data collapse is 
observed. Circles and triangles correspond to -L = 128 and 32 respectively. 
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FIG. 9: The scaling function U (s) with s = t^^^ /x is plotted on a double-log scale, for the 2D Ising 
model with the domain interface and free surface. Data collapse for different x is observed. The 
dashed line shows the power law fit, and the solid line indicates the logarithmic fit. 
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FIG. 10: The height function h{t) of the 2D Ising model is plotted on a double-log scale. The three 
upper solid lines are for the domain interface, free surface and bulk (from below). The lower solid 
lines are the pure height function Dh{t,L) of the free surface and domain interface (from below). 
Dashed lines show the power-law fits, and stars are from a double power-law fit. 
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FIG. 11: The roughness function io'^{t) of the 2D Ising model is plotted on a double-log scale. The 
three upper solid lines are for the domain interface, free surface and bulk (from above). The lower 
solid lines are the pure roughness function Duj'^{t,L) for the domain interface and free surface 
(from above). Dashed lines show the power-law fits, and stars are from a double power-law fit. 
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FIG. 12: The pure roughness function Duo'^{t, L) of the domain interface is plotted for the 2D Ising 
model with solid lines on a double- log scale, for L = 512, 256, 128 and 64 (from below). According 
to Eq. ([T9|), data collapse is observed at the curve of L' = 256. Solid circles, solid squares, and 
solid triangles correspond to L = 512, 128 and 64 respectively. 
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